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ABSTRACT 

A  model  is  presented  for  bistatic  scattering  from  ocean  sediments.  It  treats  scattering  due 
to  both  roughness  of  the  seabed  and  volume  inhomogeneities  within  the  sediment.  According¬ 
ly,  the  scattered  intensity  is  assumed  to  be  a  sum  of  two  terms,  one  proportional  to  the  rough¬ 
ness-scattering  cross  section  and  the  other  proportional  to  the  volume-scattering  cross  section. 
The  model  is  tested  against  data  acquired  as  part  of  the  Coastal  Benthic  Boundary  Layer 
(CBBL)  research  program.  As  part  of  that  program,  an  autonomous,  circularly  scanning  sonar 
system  was  deployed  in  well-characterized  regions.  This  sonar  operated  at  40  kHz,  had  a  5° 
horizontal  beam,  and  acquired  backscattering  data  over  a  50-m  radius.  During  part  of  the 
deployment,  it  operated  in  conjunction  with  a  mobile  receiving  array  so  as  to  acquire  bistatic 
data.  The  experimental  apparatus  and  procedures  are  presented,  and  results  are  compared  with 
model  predictions. 
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1.  INTRODUCTION 

The  model  for  bistatic  scattering  from  the  seabed  discussed  in  this  report  is  an  exten¬ 
sion  of  the  model  for  high-frequency  backscattering  strength  presented  by  Jackson  et  al.^ 
in  1986,  improved  by  Mourad  and  Jackson,^  and  employed  by  Jackson  and  Briggs.^  The 
model  treated  here  is  intended  for  use  at  high  frequencies  (10-100  kHz).  A  brief  discussion 
of  the  bistatic  model  has  been  given  by  Jackson,*^  and  references  to  prior  work  by  other 
investigators  can  be  found  in  Refs.  1-3. 

Backscattering  of  high-frequency  (10-100  kHz)  acoustic  energy  from  the  seabed  has 
been  the  subject  of  investigation  for  many  years.  ^  Recent  interest  in  scattering  geometries 
where  the  transmitter  and  receiver  are  not  co-located  has  led  to  the  development  of  a 
bistatic  scattering  model  that  incorporates  much  of  the  same  physics  as  the  backscattering 
model  of  Ref.  2.  In  particular,  it  includes  contributions  to  scattering  from  both  the  rough 
sediment  interface  and  inhomogeneities  within  the  sediment.  Both  the  model  and  its  incor¬ 
poration  into  an  integral  expression  for  bistatic  scattered  intensity  are  given  in  Section  2. 

Tests  of  the  bistatic  model  require  both  acoustic  experiments  in  which  specified  geom¬ 
etries  can  be  accurately  realized  and  concurrent  determination  of  the  environmental  param¬ 
eters  needed  as  inputs  to  the  model.  Such  was  the  case  in  the  experiments  described  in 
Section  3.  Both  backscattering  (monostatic)  and  bistatic  scattering  experiments  are 
described,  since  the  bistatic  model  includes  backscattering  as  a  special  case. 

The  fundamental  quantity  predicted  by  the  model  is  the  bistatic  scattering  strength. 
Section  4  compares  model  predictions  and  experimental  results  for  two  different  sediments. 
The  expression  for  the  scattered  intensity  given  in  Section  2  plays  a  fundamental  role  in  this 
comparison. 
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2.  BISTATIC  MODEL  AND  SCATTERED  INTENSITY  CALCULATIONS 

The  bistatic  scattering  strength  model  is  presented  in  Section  2.1  and  its  incorporation 
into  expressions  for  calculating  scattered  intensity  in  Section  2.2.  The  expressions  for  scat¬ 
tered  intensity  are  the  direct  link  between  the  model  and  the  experiments,  from  which  val¬ 
ues  for  bistatic  scattering  strength  are  obtained  using  a  method  defined  in  Sections  2.2  and 
4. 


The  bistatic  scattering  model  is  a  generalization  of  a  backscattering  model  that  has 
been  tested  against  data  with  good  success.  The  model  treats  scattering  due  to  both 
roughness  of  the  seabed  and  inhomogeneities  in  the  sediment  volume.  Accordingly,  the 
intensity  is  assumed  to  be  a  sum  of  two  terms,  one  proportional  to  the  roughness-scattering 
cross  section  and  the  other  to  the  volume-scattering  cross  section.  It  is  assumed  that  the 
acoustic  penetration  of  the  seabed  is  slight,  so  sediment  volume  scattering  can  be  described 
as  a  surface  process  and  quantified  by  an  effective  interface-scattering  cross  section. 

A  major  assumption  in  the  model  is  that  the  sediment  can  be  treated  as  a  lossy  fluid; 
any  effects  due  to  elasticity  or  porosity  are  neglected.  It  is  further  assumed  that  there  are  no 
gradients  in  sediment  properties,  apart  from  the  random  fluctuations  responsible  for  vol¬ 
ume  scattering.  Thus  the  sediment  can  be  characterized  by  three  parameters:  mass  density, 
sound  speed,  and  acoustic  absorption  coefficient.  The  seabed  relief  is  assumed  to  be  an  iso¬ 
tropic,  two-dimensional  Gaussian  random  process  completely  determined  by  a  spectral 
density  that  follows  a  simple  power  law  in  wavenumber.  This  adds  two  more  parameters  to 
the  model:  the  exponent  of  the  power  law  and  a  parameter  that  sets  the  overall  spectral 
level.  The  volume  scattering  strength  is  also  assumed  to  follow  a  power-law  form,  which 
adds  the  final  three  parameters  to  the  model:  the  exponent  of  the  power  law,  a  parameter 
that  sets  the  overall  spectral  level,  and  a  parameter  that  relates  density  and  compressibility 
fluctuations.  Volume  scattering  is  assumed  to  be  weak  in  the  sense  that  the  scattered  field 
is  much  smaller  in  magnitude  than  the  incident  field  (defined  as  the  field  that  would  exist 
in  the  sediment  in  the  absence  of  volume  scattering). 

The  link  between  the  bistatic  scattering  model  and  the  experiments  is  the  model  for  the 
intensity  of  scattering  from  the  seabed  as  a  function  of  time.  The  experiments  used  trans¬ 
mitters  and  receivers  with  known  beam  patterns  and  with  measured  separations,  orienta¬ 
tions,  and  heights  above  the  bottom.  The  transmitted  pulses  were  short  (2  ms).  In  the 
experiments,  the  energy  received  at  any  point  in  time  depends  not  only  on  the  bistatic  scat¬ 
tering  cross  section  of  the  seabed  but  also  on  the  geometrical  parameters  and  the  type  of 
pulses  used.  The  intensity  model  accounts  for  all  these  factors  via  an  integration  over  the 
area  of  the  seabed  that  is  scattering  energy  into  the  receiver  at  any  point  in  time.  The  results 
of  this  integration  were  directly  compared  with  the  experimental  results,  and,  where  appro¬ 
priate  (see  Section  4),  experimental  values  for  bistatic  scattering  strength  and  the  associated 
bistatic  angles  were  extracted. 
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2.1.  Bistatic  Model 

2.1.1.  Definitions 

The  equation  for  bistatic  scattering  strength  is  of  the  form 

.  (1) 

where  log  is  base  10  and  (0^,  <|)^,  Q.)  and  (6j»  6,)  are,  respectively,  the  rough¬ 

ness  and  volume  contributions  to  the  scattering  cross  section  per  unit  area.  The  angles  0^, 
<j)^,  and  0^  are  defined  in  Figure  1.  The  incident  grazing  angle  is  denoted  0^.,  the  scattered 
grazing  angle  is  denoted  0^,  and  the  bistatic  angle,  defined  as  the  difference  in  azimuth 
between  the  incident  and  scattered  directions,  is  denoted  (j)^.  In  general,  one  needs  four 
angles  to  describe  bistatic  scattering,  two  grazing  angles  and  two  azimuths;  however,  only 
the  azimuthal  difference  is  needed  here  because  bottom  statistics  are  assumed  to  be  trans¬ 
versely  isotropic.  Later  expressions  for  the  scattering  cross  section  will  employ  the  follow¬ 
ing  geometric  parameters: 

A  =  ^  [  (cos0.)^-2cos0.cos0  cos(t)  +  (cos0  (2) 

r  ^  i  i  s  ■  s  s 


and 


A  =  i  (sin0.  +  sin0  ) . 
Z  ^  s 


(3) 


z 
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These  dimensionless  parameters  are  proportional,  respectively,  to  the  transverse  and 
vertical  components  of  the  change  in  the  acoustic  wave  vector  upon  scattering.  A  parameter 
proportional  to  the  magnitude  of  the  change  is  also  used. 

A  =  jAf  +  Aj.  (4) 

The  bistatic  model  treats  the  sediment  as  a  fluid  which  is  homogeneous  except  for  the 
fluctuations  responsible  for  volume  scattering.  The  mean  values  of  fluid  properties  are 
defined  by  three  dimensionless  ratios,  p,  v,  and  8.  These  are,  respectively,  the  sediment/ 
water  density  ratio,  the  sediment/water  sound-speed  ratio,  and  the  ratio  of  imaginary  to  real 
wavenumber  in  the  sediment.  The  statistics  of  the  sediment  properties  are  assumed  to  be 
spatially  stationary.  Thus  the  ratios  defined  above  are  independent  of  position;  that  is,  there 
are  no  gradients  in  the  mean  parameters. 

The  surface  roughness  part  of  the  bistatic  model  uses  the  Kirchhoff  approximation  near 
the  specular  direction.  In  other  directions,  it  employs  the  small-roughness  perturbation 
approximation.  Bottom  roughness  is  assumed  to  be  a  spatially  stationary  random  process 
described  by  a  two-parameter,  isotropic,  two-dimensional  roughness  spectral  density: 

W{K)  =  - (5) 

(h^Ky^ 

where  ^  is  the  magnitude  of  the  two-dimensional  wave  vector  and  the  spectrum  is  normal¬ 
ized  so  that  the  integration  over  a  finite  region  in  JiT-space  gives  the  mean-square  vertical 
deviation  of  the  seabed  from  the  mean  plane  due  to  those  Fourier  components  included  in 
the  integration.  The  exponent  is  restricted  to  the  range 

2<Y2^4.  (6) 

The  parameter  ^2  in  Eq.  5  gives  the  strength  of  the  spectrum  and  has  dimensions  (length)."^ 
The  parameter  Hq  is  simply  a  reference  length  needed  to  balance  dimensions  in  the  equa¬ 
tion.  It  is  assigned  the  numerical  value  1  and  hence  does  not  appear  in  calculations. 

The  volume-scattering  portion  of  the  bistatic  model  is  an  improvement  on  the  back- 
scattering  phenomenological  approach.  The  bistatic  model  uses  perturbation  theoty  for 
volume  scattering  along  the  lines  developed  by  Ivakin  and  Lysanov,^  Hines,^  Tang,'  and 
Lyons  et  al.^  Where  the  older  backscatter  model  employed  a  single  parameter  to  quantify 
sediment  volume  scattering  strength,  the  bistatic  model  requires  three  parameters  to  char¬ 
acterize  spectra  for  inhomogeneities  in  density  and  compressibility.  The  spectrum  for  den¬ 
sity  fluctuations  is  taken  to  be  of  the  same  power-law  form  as  the  roughness  spectrum: 
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(k) 


(7) 


This  spectram  is  isotropic  and  normalized  such  that  an  integral  over  a  finite  volume  of  fe- 
space  yields  the  mean-square  density  fluctuation  divided  by  the  square  of  the  mean  density. 
Compressibility  fluctuations  are  treated  analogously  and  are  assumed  to  be  proportional  to 
the  density  fluctuations.  This  is  equivalent  to  the  assumption  used  by  other  authors.^’^  Thus 
the  spectrum  of  compressibility  fluctuations  is 

'^KK  = 


and  the  cross  spectrum  is 


'^p*r  =  ^^'^pp•  W 

where  p  is  a  dimensionless  parameter  and  is  one  of  the  eight  environmental  parameters 
needed  to  calculate  bistatic  scattering  when  using  the  present  model.  These  eight  parame¬ 
ters  are  defined  in  Table  1  and  were  described  in  the  third  paragraph  of  Section  2.  Also 
given  in  Table  1  are  the  numerical  values  used  for  these  parameters  (see  Section  3.4)  for  the 
two  sites  examined  experimentally.  With  these  definitions  in  hand,  we  are  now  ready  to  dis¬ 
cuss  the  roughness  and  volume  cross  sections  of  Eq.  1, 


Table  1.  Bottom  parameters  used  as  model  inputs. 


Symbol 

Definition 

Short  Name 

Eckemforde 

Panama 

City 

p 

Ratio  of  sediment  mass  density  to 
water  mass  density 

Density 

Ratio 

1.18 

1.97 

V 

Ratio  of  sediment  sound  speed  to 
water  sound  speed 

Sound  Speed 
Ratio 

0.991 

1.126 

5 

Ratio  of  imaginary  wavenumber  to 
real  wavenumber  for  the  sediment 

Loss 

Parameter 

0.00186 

0.0166 

73 

Exponent  of  sediment  inhomogeneity 
spectrum 

Inhomogeneity 

Exponent 

4 

4 

^3 

Strength  of  sediment  inhomogeneity 
spectrum  (cm^)  at  wavenumber 
(27c)  /a,  =  1  cm’^ 

Inhomogeneity 

Strength 

0.00013 

0.0000161 

11 

Ratio  of  compressibility  to  density 
fluctuations  in  the  sediment 

Fluctuation 

Ratio 

-0.69 

-2.44 

72 

Exponent  of  the  bottom  relief  spectrum 

Spectral 

Exponent 

3.420 

3.12 

W2 

Strength  of  bottom  relief  spectrum  (cm"^) 
at  wavenumber  (2n)/X  =  1  cm“^ 

Spectral 

Strength 

0.00231 

0.00849 
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2.1.2.  Bistatic  Roughness  Scattering 


In  this  section,  expressions  are  presented  for  the  bistatic  cross  sections  in  the  Kirchhoff 
and  perturbation  approximations.  The  net  cross  section,  0,) » appearing  in  Eq.  1 

is  formed  by  smooth  interpolation  between  the  Kirchhoff  cross  section  near  the  specular 
direction  and  the  perturbation-theory  cross  section  elsewhere.  The  interpolation  between 
these  two  approximations  is  defined  so  that  the  smaller  cross  section  takes  precedence.  This 
procedure  is  based  on  the  fact  that,  for  power-law  spectra,  the  perturbation  approximation 
overpredicts  scattering  near  the  specular  direction  owing  to  the  singularity  in  the  relief 
spectrum  at  zero  wavenumber.  In  contrast,  the  Kirchhoff  approximation  tends  to  overpre¬ 
dict  in  other  directions.^  The  interpolation  scheme  used  here  is 


(10) 


with  T|  =  -2 .  The  and  a  ^  cross-section  expressions  are  defined  next  and  are  stan¬ 
dard  results,  expressed  here  in  notation  convenient  to  the  bistatic  application. 


2. 1 .2. 1  Kirchhoff  Approximation 

Analogous  to  the  monostatic  expression  used  by  Jackson  et  al.,*  the  bistatic  Kirchhoff 
cross  section*®’^  ^  can  be  expressed  in  the  following  form: 


\R(K) 


\2r 


A2 


f 

■J 


(m)  udu, 


(11) 


where 

^  =  2k2Afc2(2kAp-2«.  (12) 

In  Eq.  11,  Jq  (m)  is  the  zeroth-order  Bessel  function  of  the  first  kind  and  k  is  the  acoustic 
wavenumber  in  water.  The  parameters  a  and  are  roughness  structure  function  parame¬ 
ters  related  to  ^2  follows: 

a  =  ^-l  (13) 

2 

and 

27tw.,r(2-a)  2-2« 

C2  =  - LJ - ^ - .  (14) 

( 1  -  a)  r  ( 1  +  a) 
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The  function  R  (Q.^)  is  the  complex  plane-wave  reflection  coefficient  (the  so-called  Ray¬ 
leigh  or  Fresnel  reflection  coefficient)  for  a  flat  interface  separating  water  and  sediment.  It 
can  be  expressed  in  terms  of  the  parameters  p ,  v,  and  6^  and  is  evaluated  at  the  grazing 
angle 

e.^  =  asin  (A) ,  (15) 

where  A  is  defined  in  Eq.  (4).  This  angle  is  not  the  actual  incident  or  scattered  grazing  angle 
relative  to  the  horizontal;  rather,  it  is  the  grazing  angle  that  would  result  if  the  seabed  were 
tipped  in  such  a  way  as  to  provide  specular  reflection  between  the  transmitter  and  receiver. 
Thorsos  (private  communication)  finds  that  this  value  gives  improved  accuracy  compared 
to  other  choices,  as  reflection  from  suitably  oriented  facets  tends  to  dominate  the  scattering 
process  near  the  specular  direction. 


2. 1 .2.2  Perturbation  Approximation 

The  bistatic  backscattering  cross  section  computed  in  the  perturbation  approxima- 
tion^“^’^^’*^  can  be  put  in  the  following  form^'^: 

=  it''|l+R(0p|2|H.fi(ep|WH'(2*:A,)-  (16) 


Equation  16  involves  the  reflection  coefficient  and  roughness  spectral  density  discussed 
earlier.  The  argument  of  the  roughness  spectral  density  is  the  “Bragg  wavenumber.”  The 
complex  function  G  is 


G  =  —  1  COS0.COS0  cosd) - +1 - . 

Vp  A  *  ^  p  J  p 


(17) 


In  Eq.  17, 


K 


1  +  j'S 
V 


(18) 


is  the  complex  wavenumber  in  the  sediment  divided  by  the  real  wavenumber  in  water  and 
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2.1.3.  Sediment  Bistatic  Volume  Scattering 


Mourad  and  Jackson^  used  a  volume-scattering  expression  similar  to  that  of  Stock¬ 
hausen^^  in  their  backscatter  model.  The  bistatic  equivalent  is  readily  obtained.  It  relates 
the  sediment  volume-scattering  cross  section,  to  the  effective  interface  bistatic-scatter- 
ing  cross  section,  (6^,  0.) ,  appearing  in  Eq.  1. 


.  |l-fig(epp|l-Hi?(9ppa, 

^bv  s’  s’  i  2itp2/m[P(ep +P(ep] 


(20) 


Perturbation  theory  is  used  to  obtain  the  volume-scattering  cross  section,  .  Adapting 
a  result  given  by  Ishimaru^^  to  the  present  situation  and  using  Eqs.  7-9  yields 


+  cos0.cose^cos(|)^  -  P  (0.)  P  (Op  pWpp  (Ak) .  (21) 

Acoustic  loss  in  the  sediment  has  been  included  by  allowing  the  wavenumber  in  the 
sediment  to  be  complex.  The  spectrum  W  is  evaluated  at  the  Bragg  wavenumber  for  vol¬ 
ume  scattering,  which  is  the  magnitude  of  me  difference  between  the  real  parts  of  the  inci¬ 
dent  and  scattered  three-dimensional  wave  vectors  (defined  in  the  sediment). 


Ak 


4Af+  {Re  [Pie.)  +P(0p]} 


1/2 


(22) 


2.2.  Scattered  Intensity 

The  experiments  described  in  Section  3  use  a  directional  receiver  and  a  directional 
source  that  emits  a  short  FM  slide.  To  determine  values  for  the  bistatic  scattering  strength 
from  these  experiments,  an  expression  for  the  bistatic  scattered  intensity  as  a  function  of 
time  was  derived  in  terms  of  the  bistatic  cross  sections  given  in  Eq.  1 .  The  expression  given 
below  is  the  scattered  intensity  lit)  of  a  short-duration,  or  “impulse,”  pulse.  This  expres¬ 
sion  must  be  convolved  with  the  signal  transmitted  (intensity  vs  time)  to  obtain  predictions 
of  the  absolute  levels  and  temporal  structure  of  the  scattered  intensity  for  any  experiment. 

The  impulse  response  is 


/«)  =  /olJ 


(0,02)2 


10 


10 


dxdy 


(23) 


(see  definitions  in  Figure  2).  Here,  a  is  the  sum  of  the  surface-  and  volume-scattering  cross 
sections  appearing  in  Eq.  1,  is  the  attenuation  of  sound  in  water  for  the  frequency  of 
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interest  in  decibels  per  unit  length,  and  b,.  are  the  beam  patterns  of  the  transmitter  and 
receiver,  respectively,  and  Dj  and  £>2  are  the  distances  shown  in  Figure  2.  The  source  is 
assumed  to  emit  a  rectangular  pulse  of  infinitesimal  length  dt  that  has  intensity  at  a  range 
of  1  m. 


z 


Figure  2.  Bistatic  geometry  and  definitions  of  variables 


The  integrations  in  Eq.  23  are  over  the  infinitesimal  area  lying  between  the  two 
ellipses,  which  have  associated  travel  times  t  and  t  +  dt  from  transmitter  to  seabed  to 
receiver.  Therefore  the  limits  on  the  x  and  y  integrations  in  Figure  2  are  actually  functions 
of  time.  Furthermore,  at  any  time  t,  the  position  x  can  be  written  as  a  function  of  y.  Thus, 
as  shown  next,  the  integral  in  Eq.  23  can  be  rewritten  as  a  pair  of  one-dimensional  integrals 
over  y. 

Using  the  variable 
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£)  =  +£>2  =  cr, 


(24) 


where  c  is  the  speed  of  sound  in  water,  the  relation  between  x  and  y  for  a  given  t  (or  equiv¬ 
alently  D)  can  be  put  in  the  form  of  an  equation  for  an  ellipse. 


£2 


(25) 


with 


a  = 


Q 


e  = 

P  =  4(jc2-Xj)  (x|-jc2  + J2_^2_£)2^  -8XjD2, 


Q  =  D 


'  [Z)2  -  (Xj  -  X^  2  _  _  j2j  2  _  4j2j2 


d2_  (jcj  -x^y 


(26a) 

(26b) 

(26c) 

(26d) 

(26e) 

(26f) 


Using  the  enlargement  of  the  (A:,y,0)  region  shown  in  Figure  2  and  Eqs.  24  and  25  allows 
Eq.  23  to  be  rewritten  as 

a  £»■ 


I{t)  =  /o(AD)  J 


-b 


dx 

(cb  b  ) 

dD 

XT 

{D.D^y 


10 


10 


dy. 


(27) 


where  dx/dD  can  be  calculated  using  Eq.  25.  Equation  27  was  numerically  implemented 
using  the  measured  beam  patterns  of  the  source  and  receiver  and  was  used  as  described  in 
Section  4  to  obtain  experimental  values  of  the  bistatic  scattering  strength.  It  is  important  to 
note  that  in  practice  two  integrals  are  calculated  using  Eq.  27,  since  x  is  double  valued  in  y; 
i.e.,  from  Eq.  25, 


X  =  X  ±a  1  . 


(28) 
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3.  BISTATIC  EXPERIMENTS 

Sections  3. 1  and  3.2  describe  the  transmitter  and  receivers  used  in  obtaining  the  back- 
scattering  and  bistatic  scattering  results  presented  in  Section  4.  Section  3.3  describes  a  typ¬ 
ical  data-acquisition  cycle.  Section  3.4  describes  the  sites  where  the  experiments  were 
conducted;  included  in  that  section  are  numerical  values  for  the  parameters  in  Table  1. 


3.1.  The  Benthic  Acoustic  Measurement  System 

The  Benthic  Acoustic  Measurement  System  (BAMS)  (Figure  3a)  operates  at  40  kHz 
and  acquires  data  on  acoustic  backscattering  from  a  circular  region  whose  radius  is  set  by 
surface  reverberation.  The  system  employs  a  planar  transmitter/receiver  array  mounted  on 
a  rotator  mechanism  at  the  top  of  a  5-m-high  tripod.  All  system  operations  are  under  control 
of  an  on-board  computer  with  a  data-storage  capacity  of  approximately  120  Mbyte.  Storage 
batteries  provide  sufficient  energy  to  operate  the  computer,  electronics,  and  rotator  for  the 
time  required  to  reach  the  limit  of  the  data-storage  capacity.  This  time  may  be  days  or 
months,  depending  on  the  rate  at  which  data  are  acquired.  In  the  experiments  reported  here, 
the  interval  between  sonar  scans  was  variously  set  at  20  min,  30  min,  and  1  hour,  and  each 


Figure  3.  Benthic  Acoustic  Measurement  System  (a)  and  mobile  receiving  arrays  (b). 
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scan  covered  a  full  360°  in  72  steps  of  5°.  The  5°  step  increment  was  chosen  on  the  basis 
of  the  horizontal  directivity  patterns  of  the  transmitter  and  receivers,  which  were  essentially 
identical  in  the  horizontal  plane  with  a  full  width  of  6.5°  at  one-half  maximum  power.  The 
full  array  aperture  was  used  for  transmission,  giving  a  vertical  directivity  pattern  with  a  full 
width  of  14°  at  one-half  maximum  power.  For  reception,  the  aperture  was  divided  into 
upper  and  lower  halves,  yielding  vertical  directivity  patterns  of  27°  for  the  receivers.  The 
levels  of  the  first  vertical  sidelobes  were  down  13-15  dB  relative  to  the  main  lobe  for  the 
various  directivity  patterns.  The  maximum  response  axis  of  the  array  was  depressed 
approximately  12.5°  below  the  horizontal.  Because  the  tower  was  not  perfectly  vertical,  the 
depression  angle  varied  slightly  as  the  transducer  was  rotated.  A  pendulum  potentiometer 
measured  the  depression  angle  at  each  step  in  each  scan.  The  azimuthal  orientation  of  the 
apparatus  was  measured  by  an  on-board  fluxgate  compass  and  by  a  magnetic  compass 
mounted  on  the  array  which  was  read  visually  by  divers  assisting  in  the  deployment.  The 
time  interval  between  rotational  steps  was  5  s,  of  which  about  2  s  were  required  to  achieve 
the  change  in  direction.  The  remaining  time  allowed  the  transmitter/receiver  array  to  come 
to  complete  rest  before  each  transmission. 

The  transmitted  signal  was  a  2-ms-long,  FM  waveform  with  a  constant  amplitude  over 
its  duration.  During  this  2-ms  interval,  the  frequency  was  swept  from  39  to  41  kHz.  The 
transmitter  provided  a  source  level  of  217  dB  re  1  |a,Pa/m.  The  transmitted  waveform  was 
digitized  and  stored  for  each  transmission.  These  data  show  that  the  transmitter  voltage 
remained  constant  to  within  a  small  fraction  of  1  dB  during  each  experiment.  Each  of  the 
two  receiver  channels  was  digitized  with  an  interval  between  samples  equal  to  one  quarter 
of  a  cycle  at  40  kHz.  Pairs  of  adjacent  samples  were  saved  every  0.5  ms  to  provide  a  pass- 
band  complex  signal.*^  This  sampling  rate  is  very  close  to  the  limit  imposed  by  the  sam¬ 
pling  theorem  for  bandlimited  signals,  and  some  aliasing  occurred  as  a  result.  Sampling 
was  started  10  ms  prior  to  signal  transmission;  these  early  samples  were  used  to  estimate 
the  noise  level.  The  sampling  window  extended  to  81  ms  after  transmission,  which  allowed 
digitization  of  approximately  1400  circular  scans  before  the  data  storage  capacity  of  the 
system  was  exceeded.  Digitization  and  signal  generation  are  controlled  by  a  single  clock. 
This  makes  it  possible  to  make  sensitive  comparisons  between  echoes  acquired  in  separate 
scans,  even  when  the  scans  are  weeks  apart.  Such  processing  has  been  used  to  obtain 
remote  temperature  measurements^^  and  to  observe  changes  in  the  seabed. The  backscat- 
tering  results  to  be  presented  in  Section  4  are  derived  from  this  data  set. 

3.2.  Mobile  Bistatic  Receiving  Arrays 

During  the  bistatic  measurements,  the  ship-mounted  mobile,  steerable  arrays  shown  in 
Figure  3(b)  were  also  deployed.  The  array  on  the  BAMS  tripod  served  as  the  transmitter, 
and  the  ship-deployed  arrays  served  as  receivers.  The  receiving  arrays  were  steered  by 
means  of  a  hydraulically  controlled  rotator,  which  can  be  seen  at  the  top  of  the  arrays  in 
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Figure  3(b).  The  hydraulic  controls  were  available  during  data  acquisition,  so  different 
geometries  could  be  obtained  within  a  data-acquisition  cycle  (see  next  section). 

The  mobile  arrays  were  divided  in  quadrants,  each  32  cm  long.  The  narrow  and  wide 
beamwidths  of  the  quadrants  were  8°  and  37°.  The  bistatic  data  were  acquired  with  the  hor¬ 
izontal  array,  with  each  quadrant  being  recorded  separately  at  two  different  gains  to 
increase  the  dynamic  range.  The  analog  signals  were  conditioned  using  the  APL  MAME 
system^®;  the  conditioned  data  were  recorded  using  Lab  VIEW  software  controlling 
National  Instruments  data-acquisition  boards.^^ 


3.3.  Data  Acquisition  Cycle 

A  simplified  diagram  of  the  bistatic  experiment  is  shown  in  Figure  4.  The  three  inde¬ 
pendent  bistatic  angles  ( 0.,  6  ,0  )  are  also  shown  in  the  figure.  The  mobile  array  was  sus¬ 
pended  from  the  tending  vessel,  which  was  placed  in  a  one-point  moor  near  the  BAMS 
tripod.  It  is  necessary  to  use  a  mobile  receiving  array  to  obtain  a  statistically  significant 
number  of  bistatic-scattering  cross-section  measurements  in  any  region  of  bistatic  scatter¬ 
ing  angle  space.  The  horizontal  distance  between  the  source  and  receivers  ranged  from  20 
to  100  m.  In  the  data  presented  here,  the  horizontal  distance  was  always  less  than  60  m. 

A  data  cycle  consisted  of  taking  data  with  the  mobile  array  during  a  360°  rotation  of 
the  tripod  array.  This  rotation  comprises  72  positions  (5°  increments)  at  each  of  which  the 
40-kHz  signal  was  emitted.  As  shown  in  Figure  4,  the  mobile  array  was  steered  so  that  the 
centers  of  the  transmit  and  receive  beams  intersected  each  other  on  the  seabed. 
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One  primary  data-analysis  effort  was  determining  the  bistatic  angles  from  the  data  set. 
Measuring  all  the  geometrical  parameters  needed  to  determine  the  bistatic  angles  obtained 
during  each  data  run  requires  the  use  of  the  split-beam  phase  from  the  receiver-array  quad¬ 
rants  (time  resolved  to  ±5  p,s)  as  well  as  several  supplemental  devices,  including  com¬ 
passes  (accurate  to  within  ±2°),  ranging  transducers  (±1  m),  an  altimeter  (±30  cm),  a 
pressure  gauge  (±30  cm),  and  inclinometers  (±0.1  ° ).  The  largest  uncertainty  in  the  exper¬ 
iment  is  in  the  ranging  since  it  could  not  be  performed  on  a  ping-by-ping  basis.  The  range 
was  acquired  an  average  of  once  every  5—10  pings.  Because  of  this,  data  taken  when  the 
receiving  array  moved  significantly  (>3  m)  between  range  measurements  were  not  pro¬ 
cessed.  The  uncertainty  in  range  is  felt  to  be  on  the  order  of  2  m  on  the  data  set  with  the 
most  receiver  movement.  The  uncertainties  quoted  for  the  altimeter  and  pressure  gauge 
come  from  side-by-side  measurements  at  10  depths  throughout  the  water  column  in  which 
the  difference  between  the  two  results  had  a  standard  deviation  of  26  cm.  The  uncertainty 
quoted  for  the  compass  is  greater  than  that  claimed  by  the  manufacturer  but  is  felt  to  be  real¬ 
istic  given  the  calibration  of  the  compass  while  embedded  in  the  receiver  apparatus  shown 
in  Figure  3.  The  inclinometers  were  also  calibrated  while  on  the  receiver  array  by  use  of  a 
digital  inclinometer  with  ±0.01  °  accuracy. 

The  bistatic  scattering  strengths  given  in  Section  4  were  derived  by  comparing  the 
pressures  received  at  the  face  of  the  mobile  arrays  with  those  predicted  by  simulations  of 
the  experiment  that  used  Eq.  27,  the  experimental  geometry,  the  transmitted  pulse  length, 
and  the  beam  patterns  of  the  receiving  and  transmitting  arrays. 

3.4.  Experimental  Sites 

3.4.1.  Eckemforde  Bay 

The  properties  of  the  sediment  at  this  site  are  documented  in  Ref.  22.  The  sediment  had 
a  very  high  water  content,  with  a  porosity  of  approximately  85%,  and  rather  low  acoustic 
attenuation  (0.18  dB/kHz/m).^^  The  most  striking  acoustic  feature  at  this  site  was  a  layer  of 
free  methane  bubbles  with  an  upper  boundary  0.5  to  2.0  m  below  the  sediment/water  inter- 
face.^^’^^  The  size  distribution  of  these  bubbles  was  measured  using  computer  tomogra¬ 
phy,^^  and  estimates  from  these  data  indicate  that  volume-scattering  strengths  should  be 
very  high,  between  -10  and  -20  dB/m.^^  Tang  et  al.^^  have  shown  that  these  gas  bubbles 
are  the  principal  cause  of  backscattering  at  40  kHz.  This  is  due  to  the  combination  of  low 
attenuation  and  high  volume-scattering  strength.  The  model  parameters  for  this  site 
(Table  1)  were  determined  as  part  of  the  CBBL  program  and,  with  one  exception  (discussed 
below),  are  those  given  in  Ref.  29.  In  situ  acoustic  probes  were  used  to  determine  the  values 
for  the  compressional-wave  velocity  and  attenuation  in  the  sediment.  Diver  cores  were  used 
to  determine  sediment  density,  and  underwater  photography  was  used  to  determine  the  sur¬ 
face-roughness  parameters.  From  comparisons  of  the  backscattering  data  with  the  back- 
scatter  model  of  Ref.  2,  it  was  determined  that  surface-roughness  effects  played  only  a 
small  role  in  backscattering  from  this  sediment  (a  result  consistent  with  Ref.  28).  The  vol- 
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ume  scattering  used  in  the  present  model  (characterized  by  Y3,  W3,  and  p)  comes  from  a  per¬ 
turbation  analysis,  and  bubble  scattering  does  not  fit  this  assumption.  However  these 
volume  parameters  can  be  regarded  as  phenomenological  rather  than  physical  and  can  be 
used  to  obtain  an  initial  estimate  of  the  bistatic  scattering.  In  this  regard,  the  W3  parameter 
quoted  in  Table  1  is  not  the  one  found  in  Ref.  29.  It  was  adjusted  to  give  good  initial  esti¬ 
mates  of  the  intensity  time  series  from  which  bistatic  scattering  strengths  are  derived  (see 
Section  4.2). 

3.4.2.  Panama  City 

This  was  a  coarse  sand  site.  The  parameters  shown  in  Table  1  are  those  given  in 
Ref.  29.  The  techniques  used  to  determine  the  environmental  parameters  were  similar  to 
those  used  at  the  Eckemforde  Bay  site.  Note  that  the  model  predictions  are  dominated  by 
the  measured  parameters  p,  v,  5,  Y2,  and  for  this  site. 
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4.  MODEL/DATA  COMPARISONS 

The  objective  of  the  model/data  comparisons  was  to  obtain  experimental  values  for  the 
bistatic  scattering  strength  and  compare  them  with  the  model  predictions.  For  each  trans¬ 
mission,  the  geometry  information  recorded  at  the  time  of  transmission,  the  attenuation  in 
the  water,  the  beam  patterns  of  the  source  and  receiver,  and  the  values  for  the  bistatic  cross 
section  predicted  using  Table  1  were  used  in  Eq.  27  to  obtain  an  impulse  response  by 
numerical  integration.  This  result  was  convolved  with  the  transmitted  waveform  to  predict 
the  amplitude  and  temporal  structure  of  the  intensity  scattered  from  the  seabed  into  the 
receiver.  This  prediction  was  then  compared  with  its  experimental  counterpart.  Two  typical 
examples  from  the  Panama  City  sand  site  are  shown  in  Figures  5a  (near  forward  scattering) 
and  5c  (near  backscattering).  Calculations  of  the  time  and  amplitude  of  the  direct  arrival 
and  the  signal  scattered  from  the  air/water  interface  are  also  shown. 


-surface 


0.02  0.04*  0.06  0.08  0.1 


0.02  0.04  0.06  0.08  0.1 

_ time  (s) _ 


Figure  5.  Model/data  comparisons  from  the  Panama  City  sand  site  and  corresponding  integrands 
of  Eq.  23.  Parts  (a)  and  (c)  show  experimental  waveforms  along  with  simulations  of  the  bottom- 
and  surface-scattered  signals  based  on  Eq.  27  (see  text  for  details).  The  vertical  axis  is  the  absolute 
level  at  the  face  of  the  receiver.  Parts  (b)  and  (d)  show  the  integrand  of  Eq.  23  for  the  geometries  of 
(a)  and  (c),  respectively.  The  transmitter  and  receiver  are  along  the  y  =  0  axis  symmetrically  placed 
around  x  =  0.  The  dark-gray  regions  indicate  integrand  values  at  least  30  dB  higher  than  the  light- 
gray  regions.  The  ellipses  indicate  regions  of  the  surface  that  contribute  4,  8,  12  ms,  etc.,  after  the 
initial  bottom-scattered  arrival. 
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The  air/water  signal  was  calculated  with  Eq.  27  using  only  the  surface-roughness-scat¬ 
tering  component  of  the  bistatic  model  and  roughness,  sound  speed,  and  density  parameters 
typical  of  the  ocean/air  surface.  This  calculation  was  intended  to  give  the  signal’s  time  of 
arrival  relative  to  that  of  the  bottom  signal  and  to  quantitatively  assess  its  possible  contam¬ 
ination  of  the  bottom  signal.  Near  forward  scattering  (Figure  5a),  the  air/water  signal 
arrives  well  after  the  bottom  signal;  near  backscattering  (Figure  5c),  the  air/water  signal  is 
of  much  lower  amplitude  than  the  bottom  signal.  Thus  contamination  of  the  bottom  signal 
is  mitigated  both  by  time  separation  and  by  beam-pattern  effects.  Similar  checks  of  other 
geometries  demonstrated  that  the  data  examined  here  were  not  significantly  affected  by  air/ 
water  interface  scattering. 

Figures  5b  and  5d  are  plots,  in  decibels,  of  the  integrand  in  Eq.  23  for  the  curves  shown 
in  Figures  5a  and  5c,  respectively.  Superimposed  are  ellipses  showing  the  region  of  the  sea¬ 
bed  contributing  to  the  signal  received  4,  8,  12  ms,  etc.,  after  arrival  of  the  initial  bottom- 
scattered  signal.  This  type  of  figure  allowed  a  quick  determination  of  whether  the  scattering 
was  sufficiently  localized  that  a  definite  set  of  bistatic  angles  (Figure  4)  could  be  associated 
with  the  time  of  the  peak  bottom-scattered  intensity.  When  this  was  possible,  the  bistatic 
angles  were  determined,  and  the  difference  between  the  theoretical  prediction  and  the  data 
was  used  to  derive  an  experimental  value  for  the  bistatic  scattering  strength.  This  method 
has  the  advantage  of  being  straightforward,  but  it  is  not  applicable  at  some  geometries.  In 
particular,  there  are  many  cases  where  the  scattering  includes  significant  contributions 
from  two  regions:  (1)  from  the  main  lobes  of  the  two  beams,  and  (2)  from  the  sidelobes  of 
the  beams  in  combination  with  the  peak  in  bistatic  cross  section  in  the  forward  direction. 
Near  backscattering,  the  region  making  a  significant  contribution  can  be  quite  large;  how¬ 
ever,  the  times  at  which  various  regions  contribute  vary  significantly,  so  time  was  used  to 
help  determine  specific  bistatic  angles  and  associated  scattering  strengths.  The  bistatic 
angles  for  which  data  were  obtained  allow  an  assessment  of  the  model’s  predictive  capa¬ 
bility  for  the  experimental  sites. 

This  method  has  an  inherent  advantage  over  sonar-equation-type  analyses  used  to 
obtain  scattering  strength.  In  a  sonar-equation  calculation,  one  must  make  some  estimate 
of  the  surface  area  that  is  contributing  to  the  signal.  Furthermore,  if  the  bistatic  scattering 
strength  is  varying  over  the  ensonified  area,  it  is  a  challenge  to  account  for  this  effect.  In 
carrying  out  the  integral  to  predict  the  intensity  time  series,  this  is  all  implicitly  taken  into 
account. 

Backscattering  is  a  special  case  of  bistatic  scattering.  Experimental  backscattering  val¬ 
ues  were  obtained  using  the  tower  as  both  a  transmitter  and  a  receiver.  In  the  following  sec¬ 
tions,  both  the  backscattered  data  and  the  bistatic  data  are  shown  for  each  site.  The  Panama 
City  site  is  examined  first,  since  the  sediment  was  more  in  keeping  with  the  model  assump¬ 
tions.  Tables  of  the  experimental  results  are  given  in  the  appendix. 
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4.1.  Panama  City  Site 

In  Figure  6,  model  predictions  of  bistatic  scattering  strength  are  plotted  as  a  function 
of  (])  and  e  for  incident  angles  of  5°,  20°,  and  60°.  All  experimental  values  for  incident 
angles  witlfin  ±2.5°  of  these  angles  are  plotted  as  points  for  comparison  with  the  model. 
Three  views  are  shown  for  each  incident  angle  so  that  data  above  and  below  the  predictions 
can  be  seen.  The  data  points  are  about  4  dB  in  diameter.  The  agreement  between  the  data 
and  the  model  predictions  is  sufficient  to  confirm  much  of  the  overall  angular  dependence 
predicted  by  the  model.  The  5°  data  show  that  the  angular  dependence  and  scattering  level 
are  well  predicted  for  small  incident  and  scattered  grazing  angles.  The  20°  data  suggest  that 
the  increase  in  scattering  predicted  near  forward  scattering  is  also  being  captured.  All  three 
data  sets  suggest  that  the  plateau  predicted  by  the  model  for  combinations  of  incident  and 
scattered  angles  that  are  not  near  forward  scattering  and  that  are  not  too  small  is  correct. 
However,  the  20°  data  indicate  that  the  plateau  predicted  by  the  model  may  be  too  high.  All 
of  these  conclusions  are  retained  when  the  full  data  set  is  examined. 

An  alternative  method  of  examining  the  data  is  shown  in  Figure  7,  in  which  all  data  for 
incident  and  scattered  angles  within  specified  ranges  have  been  grouped  and  plotted  as  a 
function  of  (}).  Model  predictions  were  then  calculated  using  0.  =  0^  for  the  lower  and 

upper  ends  of  each  range.  In  Figure  7a,  for  example,  the  upper  curve  is  for  0.  =  0^  =1*^° 

and  the  lower  curve  is  for  0.  =  0^  =  5° .  This  technique  eliminates  data  where  the  inci¬ 
dent  and  scattering  angles  are  not  similar  but  does  allow  a  simpler  view  of  the  data/model 
comparison.  This  view  of  the  data  again  indicates  an  overall  agreement  between  the  model 
and  the  data,  but  Figure  7b  still  implies  that  the  model’s  predictions  for  the  plateau  level 
may  be  slightly  high. 

4.2.  Eckernforde  Bay  Site 

The  bubbles  contained  in  the  Eckernforde  Bay  mud  have  been  identified  as  the  major 
contributors  to  the  backscattering,^®  and  it  is  physically  reasonable  to  assume  the  same  is 
true  for  bistatic  scattering.  The  model  described  in  Section  1  was  not  developed  to  address 
strong  scatterers  within  the  sediment  volume,  so  its  applicability  for  this  site  is  question¬ 
able.  It  is  still  valid  to  use  Eq.  27  as  a  starting  point  in  determining  experimental  values  for 
the  bistatic  scattering  strength.  In  this  case,  the  model  and  the  parameters  given  in  Table  1 
for  this  site  can  be  viewed  simply  as  a  complicated  means  of  obtaining  an  initial  estimate. 

Even  though  the  model  does  not  explicitly  handle  bubbles,  it  does  allow  volume  scat¬ 
terers  that  have  both  monopole  and  dipole  components. The  model  developed  exjdicitly 
for  this  type  of  sediment  structure  by  Chu  and  Tang  and  documented  by  Chu  et  al.  indi¬ 
cates  that,  for  the  frequency  being  used,  the  monopole  and  dipole  components  of  the  bub¬ 
bles  are  the  primary  contributors  to  the  scattering.  This  fact  motivated  the  plots  of 
experimental  data  and  model  predictions  shown  in  Figure  8  made  using  the  method  used 
for  Figure  7.  The  model  curves  use  the  parameters  given  in  Table  1  for  Eckernforde  Bay. 
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Figure  6.  Model  predictions  of  scattering  strength  (Eq.  1)  for  incident  angles  of  (a)  5°,  (b)  20°, 
and  (c)  60°,  shown  as  surfaces,  and  the  experimental  data  from  the  Panama  City  site,  shown  as 
points.  Three  views  are  given  for  each  plot  so  that  the  locations  of  all  data  point  can  be  seen.  See 
text  for  further  discussion. 


Though  it  is  not  useful  to  dwell  long  on  the  model/data  comparison,  examination  of  the  fig¬ 
ure  demonstrates  that,  at  the  least,  the  model  is  not  in  obvious  conflict  with  the  overall 
behavior  of  the  data. 
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(j)  <t> 


Figure  7.  Model  predictions  of  scattering  strength  compared  with  experimental  data  from  Panama 
City  site.  Model  predictions,  shown  as  curves,  are  for  the  upper  and  lower  values  of  the  grazing  angles 
given  at  the  top  of  the  figure  (see  Table  1  for  site  parameter  values). 
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Figure  8.  Same  as  Figure  7,  but  for  Eckemforde  Bay. 
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5.  CONCLUSIONS 

A  model  has  been  presented  for  describing  bistatic  scattering  from  ocean  sediments. 
The  model  has  its  genesis  in  the  backscattering  model  of  Ref.  2.  The  model  has  been  tested 
using  data  from  two  different  sites  visited  as  part  of  the  Coastal  Benthic  Boundary  Layer 
program.  One  of  the  sites  had  a  bottom  composed  of  coarse-grained  sand;  the  other  was  a 
gassy  mud.  Of  the  two  sites,  the  one  with  coarse  sand  better  follows  the  assumptions  of  the 
model.  For  this  site,  the  overall  agreement  between  the  model  and  the  experimental  data  is 
good,  with  some  indication  that  the  model  overpredicts  the  bistatic  scattering  by  a  few  deci¬ 
bels  at  some  angles  (see  Section  4.1).  A  comprehensive  list  of  the  bistatic  experimental 
results  is  included  in  the  appendix  for  both  the  sand  and  mud  sites  for  use  in  testing  alter¬ 
native  bistatic  models. 

In  the  case  of  the  gassy-mud  site,  gas  bubbles  within  the  sediment  contribute  most  of 
the  scattering.^^  The  agreement  between  the  model  and  the  data  (shown  in  Figure  8)  is  suf¬ 
ficient  for  the  model  to  be  of  practical  use  for  this  site,  even  though  in  this  case  it  should 
not  be  viewed  a  fundamental  physical  model.  The  more  physically  correct  model  of  Chu  et 
al.^^  was  compared  with  our  experimental  intensity  time  series  to  determine  bistatic  cross 
sections  referenced  back  to  the  bubble  layer  (i.e.,  refraction  and  attenuation  within  the  sed¬ 
iment  are  taken  into  account).  One  conclusion  from  this  comparison  is  that  inclusion  of 
multiple  scattering  may  be  needed  to  fully  model  the  Eckemforde  site. 
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APPENDIX 

Tables  Summarizing  Backscattering  and  Bistatic-Scattering  Results 

Table  Al.  Summary  of  Eckemforde  backscattering  and  bistatic  scattering. 
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Table  A2.  Summary  of  Panama  City  backscattering  and  bistatic  scattering. 
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